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Abstract 

In the simplest model of horizontal U(l) symmetry with one singlet added to the 
super symmetric standard model, we systematically reconstruct quark mass matri- 
ces from the low-energy data to prove that there are only two mass matrices found 
by Binetruy et. al.. The same U(l) symmetry constrains the hierarchical structure 
of L-violating couplings, from which we build radiative neutrino mass matrices 
accommodating the solar and hot dark matter neutrino masses and mixing. We 
find a few patterns of acceptable charged lepton and neutrino mass matrices, most 
of which are consistent with large tan/3 ~ mt/m& only. 



1. Introduction 

Understanding the hierarchies of fermion masses is one of the fundamental problems in 
particle physics. One way to explain the structure of fermion mass matrices is to assume 
a generation-dependent U(l) symmetry, as originally explored in [|IJ. In this approach, the 
horizontal U(l) symmetry constrains non-renormalizable couplings of quarks and leptons to 
certain SU (3) c x SU{2) L x U(l)y singlet field. The hierarchical values of Yukawa couplings are 
given by some powers of the vacuum expectation value of the singlet spontaneously breaking 
the horizontal U(l) symmetry. 



Recent investigations ||-[ I 1] have been amounted to reveal some remarkable features of 
this class of models. First of all, a potential connection to string theory was pointed out in 
]3|, ^j, |j| . In the simplest model with only one singlet field, phenomenologically allowed quark 
mass hierarchies require the U(l) to be anomalous. Then the cancellation of the anomalies 



may be the consequence of Green-Schwarz mechanism in string theories |P^| , a remarkable 
feature of which is the prediction of the mixing angle sin 2 9\y = 3/8 [13]. It was also shown 



that the horizontal U(l) symmetry may provide a solution to the Higgs mass problem (the 
/^-problem) || [7], [|. More recently, detailed analyses have been made to prove that this 



class of models indeed can account for the structure of the CKM mixing matrix || |TT|. In 
particular, for the simplest model with one singlet, only two sets of quark Yukawa matrices 
are found to be consistent with the experimental values ||11|| . In the first half of this paper, 



we will systematically reconstruct the quark Yukawa matrices starting from the experimental 
values to prove there are no other acceptable quark mass matrices. 

Contrary to the quark sector, the lepton sector is not much constrained. Our aim in the 
second half is to single out acceptable mass matrices for the charged leptons and neutrinos. 
Even if there are not yet decisive evidences for non-zero neutrino masses, some strong in- 



dications from astronomical and cosmological observations have been widely discussed [II 
Theoretical efforts to understand neutrino mass structure in the context of horizontal abelian 
symmetry were recently made by several groups [ r5| , |TE| , |TT[ | . They mostly rely on the see-saw 



mechanism to generate light neutrino masses. In this paper, we will look for R-parity violat- 
ing terms (in particular, L- violating terms) which are inherent in the minimal supersymmetric 
standard model (MSSM). Without an additional symmetry like R-parity, the MSSM allows 
the presence of L- or B- violating operators; 

ALQD , A'LLE , M'UDD . 

Proton stability imposes an extremely strong constraint on AA". On the other hand, in the 
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presence of the above L- violating operators neutrino masses can be generated radiatively. In 
general, these masses are not very much constrained due to the large freedom in the coupling 
constants A and A'. As we will see, this situation changes significantly once a horizontal U(l) 
symmetry determines the structure of the R-parity violating operator JlT], [l^] . Then, since the 
charge assignment is fixed to a large extent due to the known masses and anomaly cancellation 
the order of A, A' and the radiative neutrino mass matrices depends on a few charge values 
only. We will analyze these matrices and show that a few solutions compatible with the solar 
neutrino data and a hot dark matter neutrino mass exist. Simultaneously, the B-violating 
operators UDD can be forbidden by the U(l) symmetry to ensure proton stability. 

2. Reconstruction of quark mass matrices 

We consider the supersymmetric standard model whose Yukawa structure originates from 
horizontal abelian symmetry U(l)x- Taking the simplest form, we introduce only one singlet 
field whose [/(l)x-invariant non-renormalizable couplings to quarks and leptons determine 
hierarchical patterns of the Yukawa matrices. We assume that non-renormalizable couplings 
are suppressed by the string scale (or Planck scale) M invoking a string origin of the model. 
The charges of U(l)x can be any integer. When carries charge N, spontaneous breaking 
of U(l)x due to the vacuum expectation value ((f)) leaves as an unbroken subgroup. For 
convenience we normalize the charge of to —1 and thus the charges of quarks, leptons and 
Higgses are integers divided by N. We assume that spontaneous breaking of U(l)x occurs 
slightly below the scale M giving rise to the Cabbibo angle e = ((f)) /M ~ 0.22 in order to 



predict the experimental values of the CKM mixing angles [19( . 

The large mass of top quark leads us to assume that the (3,3)-component of the up- 
type quark Yukawa couplings is given by the renormalizable term Q3U3H2 which implies 
93 + ^3 + ^2 = 0. We use the capital letters Q, U, D, L, E and to denote MSSM superfields, 
and the corresponding small letters denote their charges under U(l)x- Depending on the 
value of tan/3, the bottom quark Yukawa term can be written as Q3D3H1 ((p/M) x where 
x = q% + c/3 + h\ is a not too large positive integer. Note that tan/3 ~ e x m t /mi } . The excess 
charge matrices = + Uj3 and yfj = q i3 + dj 3 (normalized by y^f = 0) determine the 
Yukawa matrices for the up- and down-type quarks 

where = qi — qj, etc. and the function 9 is defined as 8(x) = 1 when a; is a non-negative 
integer, and 6(x) = otherwise. As noted in || O, possible zeros in the Yukawa matrices 
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(H) are filled by non-minimal contributions to the kinetic term consistent with U(l)x- Their 
major role is to change the naive mixing angles calculated from the Yukawa matrices (|1|). In 
the canonical basis corresponding to normalized kinetic terms, the corrected Yukawa matrices 
are |§ 

Y u . r-u e \la\+\ u im\ Y u Y d . r^j e \lil\ + \d jm \ yd (c£, 

l,m l,m 

where, e.g., e' 9 " vanishes when \qu\ is non-interger. 

Taking the dominant contribution, one obtains for the down-type; Yf- = e x+n ^ where 

d J Qi3 + dj3 , when q i3 + dj3 + x is a non-negative integer 
[ Min [9 (\q.a\) + 9 (\d jm \) + 9 (q i3 + d m3 + x)] - x , otherwise. 

Here Min on the second line picks up the smallest value varying I, m from 1 to 3. The function 
9 assigns infinity when the argument is negative or fractional number and 9(x) = x otherwise. 
As observed in |TT|| , the expressions for n i3 and n 3j - are simple: 



= n% = | fe | , n% = \u j3 \ , n% = \d j3 \ . (4) 
We note that this property holds independently of the value of x assuming that the (3,3)- 



component is the largest: n 33 < riij. From the property > \yij\ 1 1| and yij = yi 3 + y 3 j one 
finds 

nij < n i3 + n 3j 

n i3 < riij + n 3j (5) 
n 3j < a i j + nr.; 

These inequalities n^f allow us to completely reconstruct the Yukawa matrices from the ex- 
perimental values of the mass eigenvalues and the CKM matrix. In refs. || pT| , two mass 
matrices are found to be acceptable. In the below we will give a systematic proof of their 
uniqueness. 

The eigenvalues of the Yukawa matrices Y determined by riy in (||) are e pi , e P2 and 1 where 

Pi = Min [nu +n 22 , n 12 + n 2 i) - p% (6) 
p 2 = Min [n u ,n 2 2, n 12 ,n 2 i\ 

in case that p\ > p 2 0, 0). The experimental values of p\, 2 extrapolated to the Planck 
scale are known to be (p", p^) = (8, 4) and (pf , pf) = (4, 2) for the up- and down-type quarks, 
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respectively ||20 | . Notice that one can fix n.22 = P2 (namely, n 22 = 4 and n 22 = 2) by reordering 
qi3, Ui3 and di 3 and by correct hierarchical values between the first and the second generations. 
Then, one gets riy > n 22 for i, j = 1, 2. Another important property is that 



Vn + 2/22 = 2/12 + 2/21 = Pi + P2 



(7) 



since the determinant of the Yukawa matrix (P should be regular and cannot change its value 
even after the kinetic correction. The regularity condition of Y{j also shows that both of two 
elements in the same raw or column of the left-upper 2x2 submatrix cannot be negative or 
fractional at the same time. From these properties and eq. (|B|), we arrive at the conclusion 
that e.g. the down-type excess charge matrices yielding n 22 = 2 is characterized by one of the 
following two patterns: 



y d n 
y d n 



913 + d 13 = 4 , y% 2 = q 23 + d 23 = 2 



2/22 



with q 23 = 0, -1, -2 



The latter condition comes from the fact that n 22 = \q2 3 \ + | c?23 j 
conclusion can be drawn also for the up-type quarks: 

5i3 + U13 = 8 , 523 + U23 = 4 ; 
1 923 1 + K3I =4 if q 2 3 + u 23 <0- 



(8) 

2 when y 22 < 0. Similar 



(9) 



The above properties are useful to reconstruct quark mass matrices compatible with the CKM 
mixing angles. 

The CKM matrix in the leading order can be parameterized in terms of three left-handed 
mixing angles as follows: 

/ 1 — s, , — Si 3 S 2 s — Si 3 + ^ 



V, 



CKM 



12 



3 12 + S 13 S 23 



qd 0- 
°12' : '23 



-s. 



13 



S23 + Sf 2 S 13 



'13 1 ^12^23 

__ q u q- 
23 °12°13 

1 



(10) 



where = Sf, — Sf, [^TJ. With the assumption < F33, an order of magnitude calculation 



can be made to express the mixing angles Sij in terms of the Yukawa elements [21]. In the 
model under consideration, the expressions become extremely simple. Iterative use of ([5]) 
reveals 



§13 — n 13 5 s 23 — n 23 

S12 = Min [rii 2 - ^22, ^11 + "-21 - 2n 22 ] 



(11) 
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where S{j ~ 



e Sij . From dJ), we immediately get s™ 3 = sj 3 = \qi3\. Comparing the CKM matrix 
(|TCi|) to the experimental data |19[ , we find two possibilities: 



S 12 ° r S 12 



s u l2 and s a 12 



s 13 ° r S 13 



s^ 3 and sf 3 



= 3, 
> 3 



s 23 ° r S 23 



"23 



or s 



d 

23 



2; 
2. 



(12) 



The above conditions can be rewritten as 



1 923 1 
1 923 1 



1*713 | = 3, 
|?13| > 3 



'12 



or s 



12 



'12 



'12 



(13) 



It is now useful to realize that the condition s\2 = 1 cannot be met by the term S12 = 
nn+n 2 i— 2n 2 2 since n n > 2n 22 (see (|7|)) and n 2X > n 22 = p 2 . Therefore, one gets n 12 = s 12 +p 2 
for S12 = 1. After some manipulation just with the down-type quark mass matrices satisfying 
(§), one finds that the second condition in (|T3|) cannot be fulfilled. In order to examine the first 
case in (|13|) , one tries out mass matrices satisfying (g i3 , q 23 ) = (±3, ±2) under the conditions 
in (H) and (|). Then we find that s" 2 or sf 2 = 1 can only be obtained for (q , i3,q , 23) = (3,2) 
and (—3,2). These corresponds to the two mass matrices found in ||11|| . The first one has 



the excess charges (gi3, #23) 
Yukawa matrix 



(3,2), (wi3,W23) = (5,2) and (d 13 ,d 



231 



[1, 0) producing the 



/ e 8 e 



(I) 



■yu 



5 £ 3 \ 



/ ,4 



e 7 e 4 e 2 



V 



Y ~ e 



1 J 



r e 3 £ 3 \ 



e 3 e 2 e 2 



V 



1 J 



For the second one, (^13, ^23) 
Yukawa matrix is 



-3,2), (ui3,u 23 ) = (11,2) and (d 13 ,d 



23 J 



(14) 

(7, 0); the resulting 



/ £ 8 e 5 e 3 \ 



(II) 



.13 



e 4 e 2 



£ U £ 2 



/ e 4 e 3 £ 3 \ 



Y ~ e 



1 / 



e 2 e 2 



V 



(15) 



1 J 



We also remark that the experimental values of the mass eigenvalues and the CKM mixing 
angles cannot be reproduced when any of j/^ is fractional. 

3. Lepton mass matrices and R-parity violating couplings 

The lepton sector is only weakly constrained compared to the quark sector. Furthermore, 
the MSSM allows L- and B-violating terms whose presence may destroy proton stability. The 
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usual wisdom is to assign R-parity to forbid both or either of them. As one can see, the 
horizontal U(l) symmetry compatible with the experimental quark and lepton mass matrices 
still has large number of free charges. As a result, it is always possible to choose some fractional 
or negative excess charges so that undesirable R-parity violating operators can be discarded. 
One can also give large enough positive excess charges to suppress the proton decay. In the 
below we will investigate the case where L-violating operators are allowed but B-violating 
operators are sufficiently suppressed. The structure of the L-violating operators as well as the 
radiatively generated neutrino masses resulting from the U(l)x symmetry will be analyzed. 
The measurement of solar neutrino deficit strongly indicates the neutrino mass and mixing; 

m„ a -2x 1(T 3 eV, 9 ea ~3x 1(T 2 ~ e 2 (16) 

which provides the right amount of resonance conversion u e —>■ v a (a = fi, r) inside the Sun 



I3|. Other hints for neutrino masses come from the deficit of atmospheric neutrinos and 
the need for hot dark matter in structure formation. Given hierarchical Yukawa matrices, 
it appears unnatural to accommodate all of three evidences which require nearly degenerate 



neutrino mass matrices \22\. Therefore we will analyze to what extent the ?7(l)-constrained 
radiative neutrino mass matrices can provide the solar neutrino mass and mixing as well as a 
hot dark matter neutrino mass m u ~ 10 eV. 

For the charged leptons, one gets the following Yukawa matrix elements: 

*3 ~e»+*& + with I;, ■ e j3 (17) 

where we have used the b-r unification at the Planck scale, x — I3 + e 3 + h\ — q 3 + d 3 + h\. 
The couplings of the R-parity violating terms are 



Ajjfc 


~ " 0(<); jh ) 


with 


Vijk 


— lo + li3 


+ Qj3 + d k3 


A' 

ijk 




with 


Vijk 


= lo + la 


+ lj3 + efc3 


A" 

ijk 


~ 9(y" jk ) 


with 


1l" 

a ijk 


= bo + u t 


3 + dj3 + dk3 



where lo = h + Q3 + d 3 = l 3 + l 3 + e 3 and 60 = ^3 + d% + d^. One has to bear in mind that 
Aftj = A'l u = due to the antisymmetric contractions. In the similar way as in the down 
quark sector, the experimental eigenvalues of the charged lepton mass matrix; pi = 4, p 2 = 2 
0] eliminate two free charges; 

e 13 = 4 - Zi 3 , e 23 = 2 - l 23 ; 
or e 23 = -2 -/ 23 with / 23 = -2,-l,0. (18) 
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From now on, we will concentrate on the first possiblitiy as it turns out that the second one 
doesnot yield proper lepton mass matrices for our purpose. The number of free charges is 
further reduced by considering the anomaly-free conditions of U(l)x- The anomalies corre- 
sponding to SU(3) 2 - U(l) x , SU(2) 2 - U(i) x , U{1)1 - U(l) x and U(1) Y - U(1) 2 X are 

A 3 = Yl( 2 <li + u i + d i) 

i 

M = E(% + y + ( h i + h 2 ) (19) 

i 

18 2 
A i = £(3® + 3^ + 3^ + l * + 2e i) + (^1 + M 

i 

2 1 72 72 1 „2\ /j2 t2\ 



A[ = £ (g ? - 2«? + d? - 1? + e?) - (Zi? - /i 



2) 



The desirable quark mass eigenvalues satisfying J2i(Qi3 + u is) — 12, J2i(Qi3 + ^3) — Sife + 
eis) = 6 is not consistent with A3 = A 2 = Ax = ||. But the MSSM with horizontal 
U(l)x symmetry may come from superstring theory which allows the Green-Schwarz anomaly 



cancellation mechanism [0 ; k 3 g 3 = k 2 g 2 = kigf oc g 2 tring - Here k{ are Kac- Moody levels of 



the gauge groups which are integers for non-abelian groups. For the anomaly cancellation, 
the ratio of the anomalies should follow A 3 : A 2 : A\ — k 3 '■ k 2 : k\. The gauge coupling 
unification in the MSSM still holds quite well even near the Planck scale, which implies 
k 3 = k 2 = 3&i/5 = 1 ]n|. Therefore, the Green-Schwarz mechanism indicates the relations 



among the anomalies: A 3 : A 2 : A\ = 1 : 1 : 5/3. 

Then the relation J2i(li3 + ^3) = Si(^3 + e i3) results in hi + h 2 = allowing the presence 
of the Higgs mass term [iH\H 2 in the superpotential. In this case, we are facing with the 
/i-problem. In the context of horizontal abelian symmetry, one may have a natural solution by 
assuming a slight deviation of the hierarchical eigenvalues like pj = 5 || or when R-symmetry 
is responsible for the hierarchy |S]]. The conditions A 3 = A 2 = 3v4i/5 together with A\ = 
leaves three charges free which are chosen to be Iq, h 3 and l 23 . The other charges are given by 

u 3 = -93 + h, d 3 = x - q 3 - hi , e 3 = x - l 3 - hi l 3 = l -x + hi, 

qs = 2 + -x-hb + l)-hqi 3 + q 23 ), hi = \x-\{b Q + l) (20) 
5 5 3 5 5 

bo = -(Z13 + ^23 + 3Z ) - 1 



considering the first case in (|T8|). Now the couplings A, A' and A" are functions of lo and k 3 . 
The radiative neutrino masses (assuming minimal soft terms) due to the couplings A and 
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A' are 



m 



/itan/3 + A, 



8n 2 m 2 



Tr [3(A 4 )MJ(A,)MJ + (A'jMj(A')Mj] , (21) 



where (Aj)jfc = Ay-fc, etc. and A is the soft-SUSY breaking parameter of the trilinear terms and 
mo is the soft mass of squarks or sleptons. Taking the leading order terms inside the trace, 
one can pick out the integers pij and p'^ as follows; 

3 87r 2 m \m T J 
with = Min[y ilm + yjkn + y? n + ytm\ ( 22 ) 

and p\j = Mm[y' ilm + y' jkn + yf n + yl m ] 

where /i ~ A ~ m — 0.1 — 1 TeV. It is understood to neglect negative or fractional y's inside 
Min. 

The expression (|22| ) is valid at the Planck scale where m^jm T ~ 1. Renormalizing Md, 
M e as well as A, A' down to the weak scale, one gets the significant effect; A/A' ~ 3 like 
rrib/m T ~ 3. Therefore the neutrino mass matrix at the weak scale is given by 



m\j ~ m[e^~ 3 + t p ^\ (23) 

where m ~ 0.4 - 0.04MeV. 

Writing m^- ~ me n « , it is simple to get the number n\- in the case where l and l + / i3 
are all positive integers. It is convenient to calculate the neutrino mass matrix in the original 
basis with non-canonical kinetic term. As Iq, Iq + /j 3 are positive integers, A i33 is present and 
gives the largest contribution to the neutrino mass matrix so that 



n u = 2/ + 



/ 2/ 13 Z 13 + Z 2 3 ^13 
^13 + ^23 2Z 23 Z 23 

V li3 J 



3 . (24) 



When m ~ 0.04 — 0.4 MeV, one can extract the following best values for Zj 3 , Iq and bo; 

*i3 = 5, / 23 = 3, / = 4(5), 6 = ^(^) (25) 

which (for fh ~ 0.4 MeV) give m„ T ~ 20(10) eV, m Vll ~ 2(1) x 10~ 3 eV and 6 eil ~ e 2 . Even 
though the mass of a hot dark matter neutrino is not precisely fixed, one hardly can allow 
one unit of deviation in Zq since it leads to two order deviation (e 2 ) in the neutrino masses 
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which is too far from the desired value. The same is also true for the solar neutrino mass. The 
resulting neutrino mass matrices are 

,13 



V 

m ~ 



V 



.15 



.13 



.10 



.10 



.11 



e e e 

We have chosen tau-neutrino as the heaviest one. In either case in 



(26) 



5D, the excess charge 



matrix and mass matrix in the canonical basis for charged leptons are given by 



/ 



V 



5\ 
3 

-1 oj 



4 
2 



/ e 4 e 4 e 



5 \ 



e 1 e 1 



(27) 



1 J 



It is interesting to note that the charged lepton mass matrix in (|27|) allows only x = 0, that 
is, large tan/3 ~ m t /mf,. This conclusion holds for either matrix (I) or (II) in (|H|) and (|I~5|). 
Furthermore, due to the fractional value of &o the B-violating operators UDD are forbidden. 
We also notice that higher dimensional operators like QQQL, uude dangerous for proton 
stability are also forbidden by the same reason. 

Now let us generalize the above calculation. Given any fractional numbers Iq, 1^, one of 
the following cases may come out. First, some of them can be fractional leading to degenerate 
matrices with a 2 x 2 nontrivial submatrix only. Matrices of this type cannot lead to the correct 
solar neutrino and hot dark matter masses simultaneously and are therefore disregarded. 
Second, A and A' are allowed to generate proper neutrino masses. Then the operators UDD 
should be completely forbidden or highly suppressed to achieve proton stability; AA" < 10~ 26 . 
A large suppression turns out to be inconsistent with the desired neutrino masses due to the last 
relation in (^). Therefore we have to choose bo fractional or smaller than — Max[-Uj + dj + d/~] 
which forbids U DD completely. As a third case one might try to allow A' and A" and forbid 
A (or make it sufficiently small) which would also ensure proton stability. As we will see, the 
wanted neutrino mass matrix does not permit this possibility. 

To examine the above cases explicitly, we carry out a systematic calculation of neutrino 
mass matrices as follows. We start by calculating neutrino mass matrices in the original basis 



as in (|22|). Then we go to the canonical basis by performing a rotation Ki for the lepton 
doublets Lj\ 



Ki 



1 

e \ll2\ 



e \h2\ 
1 



rl'ial Ahz\ 



e |iis| \ 

e l«23l 
1 



(2f 



9 



Ri 



'12 



with 



s; 



(3 



n 



i.3 ) 



'12 

1 


S 12 = " 



1 / 



'13 





1 



(l 

1 

v s" 



o \ 



Next step is to rotate the lepton doublet fields to diagonalize the charged-lepton mass matrix. 
This can be done by the rotation matrix Rf, 

I 1 -s? 9 \ / 1 

1 

Ul3 

2 



/ 



'23 



(29) 



13 



12 



as can be seen from (|TT]). Taking the leading order, the total rotation matrix R = R\Ki is 
given by Rij ~ e r ' J where r# = and 



ri2 


= r 2 i = 


= Min 


K 2 - 2, 




ri 3 


= Min 


["12- 


2 + \l 23 \ 


\hs 


r 2 3 


= Min 


["12- 


2 + \h 3 \ 


1^23 


r 3 i 


= ^13 


?"32 


= \hs\ ■ 





(30) 



In most cases of interest A gives the dominant contribution for the radiative neutrino 
mass. For case (I) in (0), one finds no more acceptable neutrino mass matrices even for some 
negative values of Iq, l^. To see this, let us first note that the resultant neutrino mass matrix 
is only two-by-two if one of Iq, Iq + 1^ is smaller than — Max[g i3 + dj 3 ] (= —4 in case (I)). 
Hence, Iq, Iq + 1^ > — Ma.x[q i3 + d^] has to be met. As one can see the smallest value (e 5 ) in 
the neutrino mass matrix results from Iq, (Iq + ^3) = —4. For this choice, the same radiative 
neutrino masses and mixing as in fl26|) can be obtained taking Z13 = 13, Z23 = 11 and Iq = —4. 
However, this leads to 60 = 3 and therefore to a too large A". On the other hand, in case (II), 
we find one more neutrino mass matrix, which has negative integers Iq and Iq + Z 2 3- When the 
tau neutrino is the heaviest, the explicit pattern is given by 

Ji3 = 19, / 23 = -3, = -4 (-3), 6 = ^) 
for m = 0.04 — 0.4 MeV. The resultant neutrino mass matrix is 



(31) 



/ ,27 



m 



,13 



V 



16 ^8 



16 \ 



.13 



.11 



(32) 



e " e e- 

As the corresponding charged lepton mass matrix in the canonical basis we obtain 



( 4 24 19 \ 

-18 2 -3 
y -15 5 j 



t £ 4 £ 24 



.19 



.18 



15 



(33) 
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Notice that this charge matrix can be consistent with x — 0, 1,2 allowing also small tan/?. 

The same calculation can be done with the muon neutrino as the heaviest. For both cases 
(I) and (II), the values (Zi 3 , Z23) = (2, —3), Iq — 7 (8) and 60 = y (y ) give rise to the following 
neutrino mass matrix in the M e -diagonal canonical basis, 

\ 



/ e 15 e 10 



V 

m ~ 



.13 



e 10 e 5 e i 



V 



.13 



.11 



(34) 



It differs from fl26|) just by an exchange of the second and the third generation. In the charged 
lepton sector this solution is characterized by 

/ 4 7 
-1 2 -3 
2 5 J 



y 



2 \ 



V 



I e 4 e 7 e 2 \ 



M 



e x 

~ e 



e 5 e 2 e 3 



e 2 e 5 



(35) 



1 ; 



For case (II), we find one more solution with (^135/23) 
yielding the neutrino mass matrix 



(22,3), l 



-7 (-6) and b = | (f) 



V 



£ 27 £ 16 



.16 



.13 



.13 



.11 



(36) 



and the excess charge and mass matrix of charged leptons 

I 4 21 22 \ 
-15 2 3 
-18 -1 j 



\ 



( e 4 



V 



.19 



.18 



.12 



1 J 



(37) 



When the muon neutrino is the heaviest, only x = is allowed as can be seen from (pop and 
(0)- 

4. Conclusion 

We examined the model of horizontal U(l) symmetry with one singlet added to MSSM, 
which explains the desired hierarchies of quarks and leptons. This model turns out to be 
quite constrained by the experimental values of quark masses and mixing. We systematically 
reconstructed the quark mass matrices to prove that there are only two acceptable patterns 
of quark Yukawa structures. 
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We extended the analysis to the charged leptons and neutrinos. Combined use of the ex- 
perimental values of the charged lepton masses and Green-Schwarz mechanism of the anomaly 
cancellation of 2/(1) reduces many free charges to three (^13/23/0) in the lepton sector. De- 
pending on their values, the presence of L- or B- violating operators can be controlled. Concen- 
trating on L- violating operators (while suppressing B-violating operators for proton stability) 
we singled out a few patterns of the charged lepton and neutrino mass matrices which ac- 
commodate the solar and hot dark matter neutrino masses and mixing. Interestingly, only 
large tan/3 = mt/m^ (x = 0) is allowed when the muon neutrino is the heaviest. When the 
tau-neutrino is the heaviest, the quark mass matrix (I) allows large tan/3 and only the pattern 
(II) is consistent with both large and small tan/? (x = 0, 1, 2). 
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